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Inclusive Higgs boson and dijet production via Double Pomeron exchange 

M. Boonekamp*, R. Peschanski^ and C. Royon* 

We evaluate Higgs boson and dijet cross-sections at the Tevatron collider via Double Pomeron 
exchange when accompanying particles in the central region are taken into account. Such inclusive 
processes, normalized to the observed dijet rate observed at run I, noticeably increase the predictions 
for tagged (anti)protons in the run II with respect to exclusive ones, with the potentiality of Higgs 
boson detection. 

1. Not only recently, Higgs boson and dijet production via Double Pomeron (DP) exchange have attracted attention. 
It has been shown fil that a non negligeable fraction of events with double colour singlet exchange features may appear. 
The calculation was based on a Pomeron model formulated as a pair of non-perturbative QCD gluons. However, the 
lack of solid QCD theoretical framework for diffraction makes a theoretical determination difficult. Predictions with 

,_^ ■ many different mechanisms appeared since then without clear consensus on the expected rate at hadron colliders. 
An evaluation of experimental possibilities j| has given the prediction of a discovery potential for the Higgs boson 
through diffraction using outgoing (anti)proton tagging and a missing mass method. However, this evaluation again 
strongly depends on the theoretical framework. 

In the present paper our aim is to give predictions based on inclusive Higgs boson and dijet production at the 
Tevatron collider via DP exchange. Indeed, so far, the predictions have been based on exclusive production with- 
out taking into account the accompanying particles in the central rapidity region. This accompanying radiation is 
unavoidable and is present in dijet production, a DP process already detected at the Tevatron ||. This is obvious 
from the observed dijet over total mass fraction spectrum which is very different from the one expected in absence of 
,_^ radiation. 

£> In fact, the evaluations contained in the present paper may have a triple phenomcnological interest, i) Taking 

into account the accompanying particles allows one to normalize the theoretical predictions to the observed dijet rate 
and thus can be translated into a more constrained prediction for the Higgs boson, ii) The inclusive cross-sections 
are expected to be larger than the exclusive production mode iii) The signal/background ratio has to be evaluated 
differently in inclusive and exclusive events. 

In the following we use the Bialas-Landshoff exclusive model for Higgs boson and heavy flavor jet production [jj as a 
starting point. We modify it in order to take into account inclusive Higgs boson and dijet production and the addition 
of light quark jets and gluon jets (not included in the analysis of Ref. M) . We show that we are able to reproduce 
up to a scale (which will be fixed from experiment) the observed distributions, in particular the dijet mass fraction 
spectrum. Using the same framework, we give predictions for inclusive Higgs boson production via DP exchange at 
Oh. the Tevatron. 

The main lesson of our study is that, if the signal/background ratio can be maintained at a reasonable level when 
associated central radiation is allowed, an interesting discovery potential for the Higgs boson particle at the Tevatron 
run II can be expected in double (anti)proton tagged experiments. It is materialized in table I for the number of 
events as a function of Mh depending on experimental cuts and decay products. 

Note that theses estimates are sizeably higher than the more recent exclusive production predictions 0. The 
inclusion of all inclusive production modes of dijets allow us to fix the estimated proportion of quarks over gluon jets. 
It is small but not negligeable, contrary to the exclusive production. Thus DP production of dijets cannot be called 
a pure "gluon factory" (cf. last paper of Ref. 12]). One reason of this difference is for exclusive production the helicity 
conserving suppression of dijet production into quark jets vs. gluon jets which is absent in inclusive production. This 
re-normalization, in turn, enhances the prediction for Higgs boson production via the it loop. 

2. Let us introduce the formulae for Higgs boson and dijet production cross-sections via DP exchange. In the 
exclusive channel O one writes; 

do%*=C B (jj T \ 8((l-x 1 )(l-x 2 )-M 2 /s) JJ i^cPvi {1-Xi) a '< exp(-2\v?)\ , 
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for a Higgs boson of mass Mjj and two heavy quark jets (of total mass Mqq), respectively. a(t) = 1 + e + a't for the 
Pomeron trajectory (e ~ .08, a' ~ ,25GeV ), x\p (> .9) are the fraction of momentum of the outgoing p and p, V\p,, 
their 2-transverse momenta, fci^, those of the outgoing quark jets, A ~ 4GeV~ 2 the slope of the Pomeron pp coupling, 
and the constants Ch, Gqq are normalizations containing various factors related to the hard matrix elements together 
with a common non-pertubative factor G 8 due to the non-pertubative gluon coupling |l]]. Thus, the ratio Ch/Cqq is 
well defined while the overall normalization is not known, and will be determined by direct comparison with data. 
For heavy quarks of mass tuq and transverse mass mri,2 the hard contribution reads 

p (1 - P ) _ im 2 Q 

F QQ(P) = ^r^5- . P = ^2- (2) 

m T1 m T2 qq 

which, up to colour factors included in the Cqq, stands for the gg — > QQ hard cross-section in the model. It is 
worthwhile to note that this cross section is proportional to ttiq, and thus is quasi zero for light quarks. This reflects 
the known zero helicity constraint of the quark jet production mechanism of Rcf. |l| (For light quarks a zero appears 
in the forward anti(proton) direction, see dijet papers in pi). 

The method we propose to evaluate the inclusive production mechanisms for Higgs boson and dijets (including light 
quarks and gluons) is simple and not restricted to the Bialas-Landshoff model, even if we use it as a starting point. The 
main driving idea is that whatever the mechanism and normalization of the diffractive part of the process could be, 
the "hard" partonic interaction in the central region depends on the probability of finding the initial partons (quarks 
and gluons) in the central region at the short time. Since the overall partonic configuration is produced initially by 
the long-range, "soft" DP interaction, we will assume that, up to a normalization, the inclusive cross-section is the 
convolution of the "hard" partons — ► Higgs boson, partons — > jets subprocesses by the probability of finding these 
partons in both Pomerons, see Fig.l. The working hypothesis is that the expected factorization breaking from soft 
color interaction will essentially affect the normalization through a renormalization of the Pomeron fluxes, which are 
not the same as in hard diffraction at HERA. 

In this framework, QCD radiative corrections is the natural origin of inclusive production H. Indeed, our ansatz re- 
markably reproduces the dijet mass fraction seen in experiment which is obviously different from exclusive production, 
see Fig. 2. 

Hence, the Higgs boson and dijet inclusive cross-sections become: 

da™ cl = G P (x<(,p)G P (x s 2 ,p) M M daT \ s _+ x f x a 2 sh 
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where xf,x| define the fraction of the Pomerons' momentum carried by the gluons involved in the hard process, 
see Fig.l, and Gp(x\ 2 ,a0, are, up to a normalization, the gluon structure function in the Pomerons extracted from 
HERA experiments, see M. p 2 is the hard scale (for simplicity kept fixed at 75GeV 2 , the highest value studied at 
HERA) . We neglected the processes initiated by quarks in the Pomeron which are negligeable in Jffl] . 

The dijet cross-section ajj depends on two "hard" cross-sections gg — > Q^'Q^ 1 ' and gg — ► gg [gf. This gives for 5 
quark flavors (i = 1 • 5) 
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27 stands for the colour factor fraction between gluon jets and quark jets partonic cross-sections ||. Note that from 
now on , the gg — ► QWqW cross section is proportional to transverse and not to rest quark masses. Thus, all 5 quark 
flavors sizeably contribute to the dijet cross-section. The quark contribution remains smaller, due to the color factor 



and to the ~ \jp% behaviour of the gluon exchange cross-section compared to the ~ l/(p^M? (i) „ (i) ) one of the quark 
exchange cross-section but is no more suppressed by the helicity constraint. 

3. The first step we perform is a comparison of our results with the measurements performed at Tevatron, in the 
CDF experiment Q . To this end, we interfaced our generator with a fast simulation of the DO and CDF detectors, 
namely SHW jL0|. We chose as gluon content of the pomeron the result of the HI "fit 1" performed in Ref. 0. Note 
that the gluon density is used only up to the normalisation of the flux. We normalize our dijet cross section prediction 
with the CDF one. 

We first compare our results for the dijet mass fraction with the measurement of the CDF collaboration W in DP 
events. A tagged antiproton with 0.035 < £p < 0.095 and \t\ < 1 GeV 2 was required where £ p>p = 1 — Xyp is the 
momentum faction of the p, p carried by the corresponding pomeron. This quantity is reconstructed using the roman 
pot detectors installed by the CDF collaboration. After the CDF cuts to tag an antiproton and the fast simulation 
of the detector, we obtain a cross section of 14.4 nb, to be compared with the CDF measurement of 43.6 nb. At 
the present stage, at least two aspects are lacking in our study : firstly, our cross-section formulae do not include 
dissociation of the outgoing (anti)protons. Secondly, some (difficult to evaluate) fraction of the proton-tagged events 
seen in experiment are in fact fake tags, coming namely from excited N* states leaving an energetic decay proton 
in the detector. Both effects contribute to enhance the measured cross-section, and thus the difference between our 
prediction and the measured value []. We thus scale up our cross-sections by a factor 43.6/14.4 ~ 3. As shown in Fig. 
2, the dijet mass fraction spectrum is well reproduced. The CDF measurement could clearly not be described without 
radiation since the obtained dijet mass fraction would peak near one, up to detector resolution effects. 

We can now give predictions for the Higgs boson production cross sections in double diffractive events, by scaling 
our results by the abovedetermined factor . The results are given in Table 1, first column. We note the high values 
of the cross-sections, which predictg more than 10 events per fb^ 1 for a Higgs boson mass below 140 GeV. 

After interfacing the generator with the fast simulation package SHW Jl0| |, we can estimate the rates which could 
be observed in the experiments. The experimental resolution and acceptances of the roman pot detectors have been 
chosen to be similar to the DO ones for dipole detectors, namely the i-resolution is 0.1^/t, ^-acceptance |i| < 0.5 
GeV 2 , ^-resolution 0.2 %, and the ^-acceptance 100% if £ > 0.04, 0% if £ < 0.01 and linear between and 100% if 
0.01 < £ < 0.04 Q. The tagging efficiency (see column 2 of Table 1) is quite good if one uses dipole detector on 
each side. To be able to trigger these events, some activity inside the central detector will be required, and we give 
in Table 1, third column, the number of events after requiring at least two jets of pt > 30 GeV []. 

To enhance the signal to background ratio, it is possible to cut on the proton and antiproton tagged energy at 
930 GeV (see Fig. 3a), on the jet topologies (the jets coming from Higgs boson events are more central) and on the 
reconstructed mass distribution using the missing mass method [pf (see Fig. 3b, and Fig. 4). We slightly modified the 
original method to partly take into account radiation and define Mh = \/£,p£,pS Jetl — 2 ^ t2 — - — - where the Ej et i 
are the leading two jets energies, E p and E p the tagged p and p energies. We notice that the missing mass method 
is not working so nicely when radiation is included. It is however still a competitive method to reduce background 
and reconstruct the Higgs boson massy. The background over signal ratio will be detailed in Rcf || . Since we obtain 
quite high cross-sections, other Higgs boson decay channels with smaller branching ratios, like H — ► t + t~ (about 
10% of Higgs boson decay, see table 1) or H — > W + W~ are of very high interest since the expected background is 
very small. This will be further discussed in Ref. fa]. 



1 Note that the CDF cross section was measured after kinematical cuts on the proton side: 0.01 < £ p < 0.03. This cut is 
difficult to reproduce using a fast simulation of the detector since it is very sensitive to the energy measurement inside the main 
CDF detector (the proton is not tagged). Thus, we did not apply this cut on our cross section. If one tries to apply these cuts 
in our fast simulation program, one obtains a cross-section two or three times smaller depending on the energy corrections. We 
chose not to apply these additional factors to be on the safe side concerning our predictions on the numbers of Higgs boson 
events. Thus, our Higgs boson event rates may be conservative by a factor 2 to 3 if one takes into account this effect. 

2 The expected luminosity is between 20 and 25 /6 _1 per experiment for run II. 

3 This allows us to enhance our signal/background ratio (the bb diffractive background without any cuts is estimated to be 
about 2.3 10 7 events per /6 _1 , and becomes about 8.3 10 4 events after those pr cuts), as shown in Fig. 3 c,d. When we compare 
Fig. 2 c,d, we also note that radiation is more important for dijet events than for Higgs boson events, since the distribution 
for dijet events is more shifted to the left after radiation than the Higgs boson one. Hence, the pr cut is more efficient after 
radiation. 

4 Due to radiation effects which escape mostly into the beam pipe, the missing mass method cannot work as it stands and is 
modified by radiation. However, it will be very useful to perform constrained fits and to have different ways of reconstructing 
the Higgs boson mass. 
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0.8 
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5.0 


160 


5.7 


0.2 


0.1 


0.0 


4.5 


170 


3.7 


0.0 


0.0 


0.0 


2.9 



TABLE I. Number of Higgs boson events for 1 fb' 1 . The first column gives the number of events at the generator level (all 
decay channels included), and the other columns after a fast simulation of the detector. The second colum gives the number 
of events decaying into bb tagged in the dipole roman pot detectors (see text), the third one requiring additionally at least two 
jets of pr > 30GeV, the fourth one gives the number of reconstructed and tagged events when the Higgs boson decays into r, 
and the fifth one when the Higgs boson decays into W + W~ (in this channel, the background is found to be negligible). 
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Figure 1 

Inclusive production schemme. 

Xi , Vi are the longitudinal and transverse 2-momenta of the 
diffracted (anti)proton, xf, the Pomeron fraction momentum 
brought by the gluons participating in the hard cross-section, ki, 
the transverse 2-momenta of the outgoing jets in the central re- 
gion from quarks, gluons or the bb decay products of the Higgs 
boson. 
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Figure 2 

Comparison of the dijet mass fraction obtained in our model 
and CDF data (open circles). 

With radiation: The shaded distribution is the dijet distribution 
after radiation and simulation of the detector. 
Without radiation: Dotted line: distribution at generator level; 
Dashed line: after simulation of the detector. 
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Figure 3 

Comparison of diffractive dijet and Higgs boson distributions. 

the Higgs boson distributions are here always on the right of 
the bb distribution. Fig (3a, up left) gives the energy distribution 
of the tagged proton or antiproton, Fig (3b, up right) the lead- 
ing dijet mass distribution, Fig (3c, down left) the mean leading 
two jet transverse energy in case radiation is added, and Fig (3d, 
down right) is the same plot without radiation. Normalisation is 
arbitrary in all figures. 
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Figure 4 

Higgs boson mass reconstruction. On top, with radiation, on 
bottom, without. The grey distribution is the result of the miss- 
ing mass method, and the dashed line is the leading dijet mass 
distribution. 



